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TECHNICAL  REPORT  SUMMARY 

The  reeeoroh  supported  by  this  grant  is  directed  towards  the  general  prob¬ 
lems  of  discrimination  and  yield  estimation.  During  the  past  year  has  evolved 
along  two  lines:  l)  analysis  of  stress  variation  within  the  crust  and  its  possible 
effects  on  "tectonic  release";  and  2)  effects  of  large  scale  mantle  structure  and 
crustal  focusing/ defocuslng  on  parameters  such  as  m*.  trace  amplitude,  period 
and  t*  over  short  distances  near  the  Fault  lass  shot  point. 

In  section  I  and  appendix  A,  the  effects  of  stress  distribution  within  the 
Great  Basin  crust  are  examined.  Detailed  studies  of  focal  mechanisms  for 
several  areas  in  the  western  Great  Basin  show  a  systematic  change  from  strike- 
slip  motion  at  shallow  depths  to  oblique  or  normal  slip  at  depth.  Explosions 
detonated  in  Nevada  show  asymmetric  Rayleigh  wave  radiation  and  significant 
Love  wave  excitation,  usually  consistent  with  a  strike-slip  mechanism.  This  has 
been  puzzling  in  an  extensional  tectonic  environment,  where  major  earthquakes 
are  known  to  have  primarily  dip-slip  motion.  However,  it  is  entirely  consistent 
with  our  observation  that  small-to-moderate  earthquakes  at  shallow  depth  in  the 
western  Great  Basin  are  predominantly  strike-slip.  An  explanation  for  this 
behavior  is  provided  in  terms  of  rotation  of  the  axis  of  maximum  compressive 
stress  due  to  increasing  overburden  pressure  with  depth. 

In  section  II  we  present  preliminary  observations  tor  teleseisms  recorded  in 
the  vleinity  of  Hot  Creek  Valley.  Broadband  digital  recordings  were  made  near 
the  FmrtLiasM  shotpoint  for  several  months  period  during  1982.  Our  preliminary 
analysis  indicates  that  amplitudes  vary  by  more  than  a  factor  of  two  over  short 
distances  (5-10  km),  however  further  analysis  is  required  before  this  can  be 
quantitatively  established.  Travel  time  anomaly  patterns  in  a  larger  region  (20- 
30  km)  suroundlng  the  Favitlass  shotpoint  are  being  examined  using  data  from 
the  U.S.  Geological  Network  operated  in  Hot  Creek  Valley.  These  show  average 
travel  time  delays  similar  to  other  areas  of  the  Great  Basin.  Relative  travel  time 


delays  in  tha  vicinity  of  the  Faultless  shotpoint  are  -0.10  relative  to  the  stations 
in  the  eastern  ranges  bounding  Hot  Creek  Valley. 
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SKCTION  I 

ESTIMATION  OF  PRINCIPAL  STRESSES  AS  A  FUNCTION  OF  DEPTH 


iTIum 


Alan  Rjrall  and  Ute  Vetter 

Surface  wave  observations  require  a  source  with  a  significant  non-isotropic 
component  for  some  large  underground  nuclear  explosions  at  both  the  US  and 
Soviet  test  sites  (Baehe,  1982).  In  early  work  on  this  subject.  Tokaoz  and  Kehrer 
(1972a.  b)  studied  Rayleigh  wave  radiation  patterns  and  Love/Rayleigh  ampli¬ 
tude  ratios  for  a  number  of  NTS  and  Amchitka  explosions,  and  ooncluded  that 
the  non-isotropic  or  "tectonic  release"  component  was  consistent  with  a  strike- 
slip  mechanism.  For  the  USSR,  surface  waves  from  some  of  the  eastern  Kazakh 
explosions  are  more  strongly  perturbed  than  has  been  observed  for  the  US  shots 
(Baohe,  1982).  Contrary  to  the  NTS  model,  these  observations  are  consistent 
with  a  source  mechanism  aomposed  of  an  explosion  plus  a  thrust-fault  double 
couple.  For  this  source  Rayleigh-wave  phase  reversals  occur  in  all  azimuths 
when  the  seismic  moment  of  the  double  couple  is  more  than  twice  the  moment 
of  the  shot,  and  in  two  quadrants  when  the  ratio  of  the  moments  is  0.5-2.0. 

Arohambeau  (1981)  concludes  that  non-isotropic  components  would  be 
expected  to  perturb  explosion  signals  whenever  the  shot  is  detonated  in  a  pres- 
tressed  medium  or  in  the  vicinity  of  stress-loaded  faults.  He  suggests  that  the 
current  state  of  knowledge  on  tectonic  release  effects  could  be  considerably 
expanded  by  a  number  of  research  initiatives,  including  research  aimed  at 
predicting  surface-wave  perturbations  for  explosions  detonated  in  "realistic” 
prestress  environments. 

In  another  area  of  the  TBT  verification  research  program.  Archambeau 


(1962)  points  out  that  the  capability  to  identify  seismic  events  in  terms  of  their 
physical  characteristics  is  fundamental  to  the  discrimination  of  earthquakes 
and  explosions.  "In  particular,"  he  states,  "one  may  use  the  location  of  the 
•vent  along  with  knowledge  of  the  tectonics  for  the  area  in  which  the  event  has 
occurred,  obtained  in  part  from  previously  observed  earthquakes  and  in  part 
from  other  geophysical  information  including  plate  tectonics,  to  infer  the  most 
probable  (or  possible)  fault  mechanism  and  hypocentral  depth  for  an  event.” 

As  part  of  our  program  of  research  on  earthquakes  and  explosions  in  the 
western  Great  Basin,  we  have  made  a  detailed  investigation  of  changes  in  focal 
mechanism  as  a  function  of  focal  depth.  The  results  of  this  study  indicate  that 
earthquake  mechanisms  in  this  region  are  predictable  in  terms  of  epicentral 
location  and  depth.  Maximum  and  minimum  principal  stresses  can  also  be 
estimated  as  a  function  of  depth.  One  significant  result  of  the  analysis  is  that 
strike-slip  mechanisms  for  shallow  earthquakes  —  and  tectonic  release  —  at  the 
Nevada  Test  Site  are  consistent  with  the  known  extenslonal  tectonic  regime  of 
the  region.  Extension  of  this  analysis  may  provide  useful  information  on  pres- 
tress  environments  in  regions  where  other  test  sites  are  located. 

In  a  study  of  earthquakes  in  central  Nevada  (Ryall  and  Vetter,  1982)  we 
observed  that  events  with  strike-  or  oblique-slip  mechanisms  had  average  depth 
of  about  ?  km,  wheras  events  with  normal  slip  had  mean  depth  of  about  12  km. 
Subsequent  analysis  of  more  than  130  earthquakes  distributed  over  most  of  the 
western  Great  Basin  show  a  consistent  pattern  of  primarily  strike-slip  motion  for 
•vents  shallower  than  6  km,  while  those  deeper  than  9  km  generally  have  a 
strong  normal-slip  component.  Figures  1  and  2  illustrate  this  effect  for  the 
Mammoth  Lakes  area 

In  a  paper  acoepted  tor  publication  in  the  Journal  of  Gtophysical  Research. 
(Appendix  A;  Vetter  and  Ryall,  19B3)  we  explain  the  observed  change  in  mechan- 
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ism  by  increasing  overburden  pressure  with  depth.  Shallower  than  8  km  the 
maximum  compressive  stress  is  horizontal,  leading  to  strike-slip  motion. 
Deeper  than  9-10  km  the  maximum  stres  is  the  overburden  pressure,  which  is 
vertical.  At  depths  of  8-10  km  where  the  transition  takes  place  the  difference 
between  greatest  and  Intermediate  principal  stress  must  be  small,  as  indicated 
by  the  occurrence  of  both  strike-  and  oblique-slip  events  in  the  8-9  km  range 
and  oblique  to  normal  slip  at  larger  depth.  Our  observations,  together  with 
measured  stress  in  crustal  rocks  (Sibson,  1974;  Byerlee,  1977;  Zoback  and 
Zobaok.  1980;  McGarr,  1980)  can  be  used  to  roughly  estimate  the  principal 
stresses  at  greater  depth,  as  illustrated  by  Figure  3.  Within  the  western  Great 
Basin  we  find  the  maximum  and  minimum  compressive  stress  at  depth  of  3  km 
to  be  about  104  and  55  MPa,  respectively.  At  depth  of  80  km  the  maximum 
stress  is  vertical  and  equal  to  the  overburden  pressure,  about  530  MPa;  the 
minimum  stress  is  horizontal  and  about  880-300  MPa. 

Our  observations  also  indicate  that  the  fault-plane  dips  decrease  with 
depth,  but  the  change  from  strike-  to  dip-slip  motion  is  not  consistent  with  lis- 
tric  faulting  such  as  observed  in  the  eastern  Great  Basin.  Instead,  based  on  a 
model  proposed  by  Hill  (1977)  to  explain  the  predominance  of  strike-slip 
mechanisms  in  geothermal  areas,  we  suggest  that  lithospheric  extension  in  the 
western  Great  Basin  may  occur  through  oblique  or  normal  faulting  at  mid- 
crustal  depths  and  the  formation  of  vertical  fissures  or  magma-filled  dikes  at 
shallow  depth.  A  model  for  the  shallow  earthquakes  in  the  Mammoth  Lakes  area 
is  shown  on  Figure  4. 
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Figure  Captions 


Figure  1.  Locations  and  fault-plane  solutions  for  Mammoth  Lakes  earthquakes 
with  depth  less  than  9  km. 

Figure  2.  Locations  and  fault-plane  solutions  for  Mammoth  Lakes  earthquakes 
with  depth  greater  than  9  km. 

Figure  3.  Estimated  principal  stresses  as  a  function  of  depth.  Pq  —  hydrostatic 
pore  pressure,  Sj  -  —  twice  the  maximum  shear  stress  measured  for  hard 

rooks  and  extrapolated  downwards,  S*=Sa  —  smallest  horizontal  principal 
stress,  Sy  -  overburden  pressure,  Sg  —  greatest  horizontal  principal  stress. 


Shaded  area  shoes  range  of  5|.  Points  with  numbers  are  described  in  the  text. 

Figure  4,  Model  to  explain  strike-slip  faulting  for  shallow  Mammoth  Lakes  earth¬ 
quakes.  Top  left  •  Rose  diagram  showing  orientation  of  T-axis;  tope  right  -  Rose 
diagram  showing  stike  of  the  two  planes  of  the  fault-plane  solutions:  bottom  - 
model  with  vertical  Assures  (dikes)  forming  together  with  conjugate  right-  and 
left-lateral  shears  on  planes  shown. 


SECTION  II 


Teieseiamic  obtryitiona  in  the  vicinity  at  the 
Faultless  test  site 


Keith  Priestley 


Introduction 

With  the  current  emphasis  on  TTBT,  improving  the  accuracy  of  yield  deter* 
mination  to  verify  treaty  compliance  haa  given  new  importance  to  the  problems 
of  magnitude  bias  and  regional  differences  in  attenuation.  As  pointed  out  by 
Baohe  ft  of, (1981)  and  Bache(1982),  while  it  is  generally  accepted  that  there 
exist  broad  regional  (l.e.,  WUS,  EUS)  differences  in  attenuation,  there  are  still 
obstacles  to  accurate  prediction  of  amplitude  variations  at  individual  stations 
within  these  regions.  In  a  study  of  amplitude  and  m*  bias  at  LASA,  Chang  and 
von  Seggern  (1980)  found  fluctuations  of  about  0.15  m*  in  standard  deviation 
among  the  subarrays,  even  after  determining  azimuthal  corrections.  Detailed 
modeling  could  not  identify  the  source  of  the  bias  (azimuthal  effect,  subarray 
effect  or  azimuth-subarray  interaction),  and  the  authors  concluded  that  both 
amplitude  and  traveltime  anomalies  might  be  due  to  local  focusing/ defocusing 
effeota  in  the  crust.  The  authors  concluded  that  "in  a  region  as  complex  as 
LASA",  —  a  region  originally  selected  for  its  geologic  uniformity  (Green  ef  of., 
1980),  —  a  factor  of  2  in  uncertainty  results  for  predictions  at  a  station's  ampli¬ 
tude  to  a  hypothetical  event  only  50  km  away  from  a  hypothetical  calibration 
shot. 

In  a  recent  SDCS  experiment,  Teledyne-Geotech  installed  and  operated 
high-quality  short-period  stations  at  WUS  test  sites,  and  compared  the  recorded 


teleseismic  signals  with  data  from  RXON  and  two  stations  in  Maine  (Der  at  o l., 
1981).  Scatter  in  the  data  was  large  for  station  pairs  with  10-20  km  spacing.  At 
NTS  the  standard  deviation  of  relative  magnitude  differentials  was  0.2-Q.3  m*. 
However,  with  appropriate  corrections  for  crustal  structure,  the  NTS 
differentials  could  be  reduced  to  a  few  hundredths  of  a  ms  unit.  After  correc¬ 
tion,  m*  differentials  between  Great  Basin  test  sites  (NTS,  FauUlaaa  site.  Shoal 
site)  and  RKON  still  ranged  from  0.16  to  0.44  m*.  Such  large  differences  were 
not  observed  for  t*  differentials  computed  for  pairs  of  stations  by  spectral 
ratios— the  spread  for  12  WUS  stations  was  only  about  ±0.05,  while  the  difference 
between  RKON  and  the  WUS  t*  values  was  0.20.  In  agreement  with  the  study  by 
Cahng  and  von  Seggem,  Oer  at  at  suggested  that  t*  differentials  reflect  true 
regional  differences  in  anelastic  attenuation,  while  the  large  scatter  in  and 
trace  amplitude  measurements  within  the  WUS  may  be  the  result  of  local  focus¬ 
ing  effects. 

One  of  the  sites  of  interest  Is  Hot  Creek  Valley  (CNTA),  in  central  Nevada, 
the  location  of  the  FauUlaaa  test  of  January  19,  1968.  The  ISC  gave  the  FauUlass 
explosion  a  teleseismic  body-wave  magnitude  of  6.3  based  an  44  readings.  This 
is  more  than  0.3  m*  units  greater  than  the  estimated  maximum  magnitude 
(Evernden.  1971)  of  the  device  baaed  on  the  announced  yield  range  of  200  kt  -  1 
Mt  (Springer  and  Kbinamann.  1971)  Indicating  that  CNTA  has  low  attenuation 
relative  to  NTS.  However  the  Teledyne-Geotech  study  (Der  at  at,  1961)  of  SDCS 
teleseismic  data  indicates  that  relative  to  NTS,  CNTA  is  a  high  attenuation  site, 
thus  presenting  an  apparent  violation  of  reciprocity. 

Coo  logic  Betting 

Hot  Creek  Valley  in  central  Nevada  is  a  region  of  extensive,  late  Cenosoic 
voieanism.  The  voioanio  deposits  of  the  Lunar  Craters  area  in  southern  Hot 
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Creek  Valley  have  bean  studied  by  Scott  and  Ttask  (1971).  The  volcanic  field 
Includes  cinder  cones,  maars,  and  basalt  flows  of  probable  Quaternary  age. 

Structures  in  this  region  are  relatively  simple  and  are  reflected  by  the  phy¬ 
siography  typical  of  the  Basin  and  Range  province  (figure  1).  Folding,  where 
present,  is  mild,  and  numerous  normal  faults  separate  blocks  which  are  tilted 
generally  eastward.  In  plaees,  this  pioture  is  made  more  oomplex  by  large 
masses  of  Paleoxoic  rocks  that  have  slid  or  been  thrust  over  the  Tertiary  ignim- 
britas. 

The  Quaternary  basalt  flows,  cinder  cones,  and  craters  are  underlain  and 
partly  surrounded  by  Tertiary  rhyolite  tuffs  and  flows,  and  Paleoxoic  carbonate 
rooks.  Basalt  flows  and  cinder  cones  of  Miocene  age  are  the  most  widespread  of 
all  the  deposits  and  extend  from  the  northern  to  the  southern  edge  of  the  vol¬ 
canic  field.  There  is  little  infilling  of  the  moat  recent  crater,  and  the  flanks  of 
the  cone  are  not  gullied  by  erosion.  The  crater  rims  have  remained  relatively 
sharp,  and  loose  ejecta  form  a  slope  of  39* ,  a  maximum  for  cones  in  this  area, 
indicating  a  very  young  age  (at  most  several  thousand  years). 

Basalt  flows,  cinder  cone  chains,  mounds,  craters,  and  fissures  rimmed  with 
pyroclastic  ejecta  were  produced  along  normal  fruits  and  fissures  resulting  from 
regional  extension  that  started  about  17  m.y  ago,  and  has  continued  to  the 
present.  At  least  three  episodes  of  late  volcanism  occurred  throughout  the  vol¬ 
canic  field,  separated  by  brief  periods  of  relative  dormancy.  These  interludes 
may  not  have  been  distinct  and  probably  represent  overlapping  pulses  of 
activity.  Bach  episode  began  with  the  eruption  of  pyroclastic  material  that 
termed  small  mounds;  as  these  grew  larger  and  became  cons  shaped,  the  initial 
extrusion  of  highly  vesicular  basalt  flows  and  ejeota  produced  composite  or 
mixed  volcanic  structures.  This  early  explosive  phase  probably  resulted  from  a 
combination  of  the  high  volatile  content  of  the  magma  and  greater  amount  of 


ground  water  available  at  the  onset  of  volcanism.  In  the  final  stages,  as  magma 
continued  to  rise  in  the  vents,  many  cones  were  breached  where  the  lavas 
overflowed  or  eroded  the  walls  of  cinders  and  ashes.  The  weakness  of  the 
rtngwalls  was  apparently  conditioned  by  linear  fissures  and  faults  along  which 
the  oones  developed,  as  shown  by  the  preferred  northeast-southwest  orientation 
of  the  ruptures. 

Digital  recordings 

To  aid  in  explaining  the  apparent  discrepency  in  the  amplitude  of  outgoing 
vs.  incoming  waves,  the  University  of  Nevada  moved  a  digital  array  previously 
operating  in  the  Catsldll  Mountains,  New  York,  (Pomeroy,  I960)  to  Hot  Creek  Val¬ 
ley  in  mid-January,  1982.  The  locations  of  the  stations  are  shown  in  figure  1. 
The  site  labled  GROUND  ZERO  is  at  the  faultless  shotpoint;  the  station  labled 
MOREY  MINE  is  in  competent  rhyolitic  rocka  in  the  range  S  1cm  to  the  west  of 
GROUND  ZERO ;  the  station  labled  NEEDLES  is  in  volcanic  flows  10  km  east  of 
GROUND  ZERO.  The  stations  were  installed  to  study  variations  in  such  parame¬ 
ters  as  m* ,  trace  amplitude,  period  and  t*  over  short  distances  within  the  vicin¬ 
ity  of  the  faultless  shot  point,  and  to  determine  to  what  extent  these  variations 
can  be  explained  by  large  structures  within  the  mantle  (e.g.,  high-velocity  zones 
such  as  proposed  under  NTS  calderas  by  Spence,  1974),  or  by  local 
foeusing/defocusing  within  the  crust  (e.g.,  such  as  the  crustal  inhomogeneities 
at  LASA  proposed  by  Chang  and  von  Seggern,  1980). 

Several  unanticipated  problems  developed  which  limited  the  scope  of  this 
experiment.  Meet  critical  of  these  problems  involved  transmitting  the  seismic 
data  from  the  remote  stations  to  the  central  recording  site.  As  the  array  was 
operated  In  New  York,  data  was  transmitted  from  stations  to  recording  site  via 
leased  telephone  Ones,  something  not  possible  in  remote  areas  at  Nevada.  Thus 


the  system  was  modified  to  employ  VHT  radios.  The  Seismologies!  Laboratory 
technical  staff  have  had  a  great  deal  of  experience  with  VHT  transmission  of 
selsmio  data  over  the  past  eight  years  and  anticipated  no  problems  with  the 
relatively  short  transmission  paths  involved  in  Hot  Creek  Valley.  However,  we 
experienced  a  great  deal  of  intermittent  radio  interference  especially  in  the 
path  from  ground  zero  to  the  central  recording  site.  At  times,  three  watt  radios 
were  insufficient  over  short  (15*80  km)  line-of-site  paths,  while  the  technical 
personnel  were  simultaneously  maintaining  visual  contact.  We  have  not  experi¬ 
enced  similar  interference  at  any  other  sites  in  western  Nevada  and  suspect  that 
the  interference  may  result  from  the  Air  Force  installations  in  the  valley  and 
surrounding  ranges.  Thus,  over  much  of  the  seven  months  of  operation,  data 
were  received  from  tin  remote  sites  only  intermittently. 

Other  problems  included  poor  environmental  controls  for  the  recording 
computer  facility  (heat,  dust,  power  fluctuations,  lightning),  and  Inexperienced 
local  peraonneL 

Much  time  during  the  past  year  has  involved  writing  and  assembling 
software  for  playing  back  and  processing  the  digital  data.  We  have  made  a  prel¬ 
iminary  analysis  of  the  two  seismograms  shown  in  Figure  2.  In  both  instances 
the  amplitude  observed  at  the  "ground  zero"  site  are  significantly  larger  than 
thoee  observed  at  the  two  stations  sited  on  rock.  All  stations  were  operating  at 
the  same  magnification  and  hence  the  amplitude  variations  shown  in  the  time 
series  between  components  for  an  event  represents  real  variations  in  amplitude 
between  the  three  sites.  No  corrections  have  been  made  for  site  amplification 
resulting  from  differences  in  sediment  thickness  beneath  each  site. 

Beth  events  were  from  the  northwest,  and  both,  qualitatively  stow  a  high 
degree  of  phase  ooherenoe  over  the  first  several  cycles  of  the  P-wave  on  all  com¬ 
ponents.  The  first  event  Is  larger  resulting  in  a  higher  signal  to  noise  ratio.  All 
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three  north-south  component*  for  this  arrant  show  a  strong  phase  following  the  P 
onset  by  5.8  ±02  seconds.  This  may  be  a  P  to  S  conversion  at  a  boundary 
beneath  the  array.  The  average  P  and  S  wave  velocities  for  the  Great  Basin 
orust  are  6.27  and  3.65  k/s  respectively  (Hill  and  Paldser,  1967;  Priestley  and 
Brune,  1978).  If  the  observed  phase  is  assumed  to  be  a  P  to  S  conversion  at  the 
base  of  the  crust,  the  boundary  occurs  at  a  depth  of  about  49  km.  This  is  deeper 
than  the  Moho  depth  postulated  within  east-central  Nevada  from  refraction  (37 
km.  Eaton.  1963;  38  tow,  Hill  and  Paldser,  1967;  35  km,  Priestley  at  al,  1982). 
The  average  P  and  S  wave  velocities  for  the  entire  lithosphere  are  6.96  and  4.04 
k/e.  respectively  (Burdick  and  Helmbarger,  1978;  Priestley  and  Brune,  1978).  If 
the  conversion  is  assumed  to  occur  at  the  base  of  the  lithosphere  the  boundary 
occurs  at  54  km  depth,  less  than  the  average  depth  for  the  base  of  the  litho¬ 
sphere  from  published  body  wave  (65  km,  Burdick  and  Helmberger,  1978)  and 
surface  wave  (64  km,  Priestley  and  Brune,  1978)  studies.  This  postulated  con¬ 
verted  phase  is  apparent  on  only  one  of  the  two  events  examined  to  date. 

figure  3  shows  a  comparison  of  the  power  spectra  of  the  vertical  com¬ 
ponent  seismograms  for  the  two  events.  The  sample  rate  of  the  data  is  25 
samplas/seoond.  and  the  time  window  taken  is  the  first  10.24  seconds  of  the  P- 
wsive.  To  estimate  the  mdse  level,  10.24  seconds  of  the  time  series  preceding 
the  P-onset  were  analysed  in  the  same  manner  as  the  P-wara.  For  both  events, 
the  amplitude  of  the  signal  at  1  Hi  observed  at  station  "Ground  Zero"  is 
significantly  larger  than  that  observed  at  the  other  sites.  Instrument  and  site 
corrections  need  to  be  made  before  a  more  quantitative  statement  can  be  made 
ooooerning  local  offsets  on  m*  estimates. 


Ttiiwlmk  P-deUys  at  Hot  Crook  Valia/ 

Evernden  (1977)  summarized  past  studies  of  station  magnitude  bias,  and 
found  that  magnitude  bias  can  be  correlated  with  travel  time  anomalies,  crustal 
thickness  tinder  the  station,  a  low-velocity  channel  in  the  upper  mantle,  high 
hoot  flow,  and  upper  mantle  velocity  structure.  Chang  and  von  Seggern  (1980) 
found  that  amplitude  anomalies  at  LASA  are  linearly  related  to  the  amount  of 
travel  time  anomaly,  implying  that  both  effects  are  due  to  a  similar  cause,  which 
they  suggest  to  be  crustal  focusing.  With  this  in  mind,  we  have  begun  an  analysis 
of  teleseiamic  P-wave  delays  in  the  vicinity  of  the  Favitla ss  shotpoint  in  Hot 
Creek  Valley.  From  late  September,  1969  through  September,  1970,  the  U.S. 
Geological  Survey  operated  seven  short-period  seismograph  stations  in  and 
around  Hot  Creek  Valley.  Each  station  was  equipped  with  a  vertical  component, 
short-period  (1  Ha)  seismometer.  The  signals  were  telemetered  to  a  central 
recording  site  and  recorded  along  with  WWVB  time  code  on  16  mm  develocorder 
film,  fflth  the  common  time  base,  dock  corrections  are  unnecessary.  The  sta¬ 
tion  locations  are  shown  in  figure  1. 

Teleseiamic  P-wmve  arrival  times  have  been  read  directly  from  the  develo- 
o order  Aims  for  41  events  in  the  distance  range  95*  to  100*.  Most  events  were 
from  the  southwest  Pacific,  the  J apan- Aleutian  arc,  or  Central  and  South  Amer¬ 
ica;  however,  a  number  of  events  were  Included  from  Europe  and  the  mid- 
Atlantic  rise.  The  quality  of  the  data  varied  significantly  depending  on  the 
amplitude  of  the  arrivals  and  the  level  of  instrumental  noise.  In  general,  the 
waveforms  showed  a  high  degree  of  coherency  across  the  network.  Often  one  or 
more  stations  were  Inoperative  and  data  were  only  included  In  the  analysis  for 
which  recordings  were  available  for  at  least  tour  stations. 

A  total  of  872  arrival  times  have  bean  read  of  coherent  phases  across  the 
array.  Occasionally  the  first  arrival  time  was  read,  but  more  often  the  time  of 
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the  first  peak  or  trough  was  read.  Readings  were  made  to  a  precision  of  0.01 
second.  To  account  for  the  uncertainties  due  to  noise,  the  phase  times  were 
bracketed  by  the  time  for  which  we  could  determine  with  a  high  degree  of  cer- 
tainity  that  the  phase  had  not  arrived,  and  the  time  for  which  we  were  certain 
that  the  phase  had  arrived.  The  error  bounds  varied  from  ±0.02  seconds  for 
impulsive  events  during  quiet  intervals,  to  ±0.15  seconds  for  weaker  arrivals 
during  more  noisy  intervals.  Typically  the  error  bounds  were  approximately 
±0.05  seconds. 

Elevation  corrections  were  applied  to  the  data  correcting  the  arrival  times 
to  those  which  would  have  been  recorded  at  the  Fault Itss  shot  point  elevation. 
In  general  this  elevation  change  was  small  (i.e.  the  maximum  elevation  change 
was  1328m).  thus  significant  errors  should  not  be  introduced  into  the  data  due 
to  the  wrong  choice  of  velocity  for  the  surficial  layer.  Absolute  P  wave  residuals 
were  calculated  with  respect  to  the  Herrin  travel-time  tables  (Herrin.  1988) 
using  the  event  location  and  origin  times  from  the  Preliminary  Determination  of 
Epicenter  bulletins.  Relative  residuals  were  then  recorded  for  each  station  with 
respect  to  the  weighted  mean  travel  time  residual,  Le. 

where  Hq  is  the  relative  residual  for  the  i*  event  at  the  j**  station,  Ran  is  the 
residual  compared  to  the  theoretical  Herrin  travel  times,  N  is  the  number  of  sta¬ 
tions  recording  the  Xth  event,  and  Ify  is  a  weighting  factor  for  the  event  at 
the  jm  station.  Not  all  stations  recorded  every  event,  therefore  the  weighted 
mean  travel  time  may  be  computed  by  a  different  subset  of  stations  for  various 
events. 

The  average  travel  Urns  difference  across  the  seven  stations  is  greater  than 
2  seconds,  typical  for  the  Great  Batin.  The  pattern  of  individual  station  retidu- 


ala  relative  to  this  average  roughly  correlates  with  the  geologic  patterns  in  the 
vicinity  of  the  stations.  The  average  station  residuals  in  seconds  are  shown  on 
Figure  1.  Station  CNBC  is  located  in  the  valley  center  and  shows  the  largest 
positive  residual  (0.23  sec),  presumably  resulting  from  the  slower  velocity  valley 
sediments.  The  stations  along  the  eastern  ranges  (CNEC,  CNHM,  CNSB,  and 
CNPS)  all  show  small  residuals  {-0.02±0.02  sec).  The  stations  along  the  western 
range  (CNHC  and  CNHR)  show  larger  negative  residual  (-0.10  sec). 

Figure  4  shows  the  distribution  of  relative  residuals  with  azimuth  and  on  an 
equal-area  lower-hemisphere  stereographic  projection.  Station  CNBC  has  all 
positive  residuals  and  shows  no  obvious  azimuthal  trend.  Stations  CNHM  and 
CNEC  generally  show  positive  relative  residuals  for  waves  approaching  at  steep 
angles  from  the  southeast,  and  negative  residuals  for  waves  approaching  from 
the  southwest.  Station  CNHR  and  CNHC,  both  located  on  the  west  side  of  the  val¬ 
ley,  have  negative  residuals  (with  a  few  exceptions  to  the  southeast),  with  no 
apparent  azimuthal  trend.  Station  CNSB  has  mainly  positive  residuals  for  most 
azimuths,  with  the  largest  residuals  for  the  southern  quadrant.  Station  CNPS 
shows  the  strongest  azimuthal  trend,  with  positive  residuals  for  waves  approach¬ 
ing  from  the  north  or  south,  and  negative  residuals  for  waves  approaching  from 
the  east  or  west.  For  all  stations,  waves  approaching  from  the  southeast,  the 
azimuth  of  the  recent  volcanism,  show  the  largest  relative  delays,  although  each 
station  has  its  own  azimuthal  pattern.  We  are  currently  timing  the  last  six 
months  of  data  from  the  U.S.  Geological  Survey  network  and  expect  that  the 
remaining  data  will  help  clarify  the  underlying  cause  of  the  delay  pattern. 

To  supplement  the  teleseismic  recordings  in  Hot  Creek  Valley,  teleselsmic 
travel  time  residuals  have  also  been  computed  from  recordings  of  the  Fauttlwss 
explosion.  Arrival  times  were  taken  from  the  ISC  bulletin,  using  only  those 
arrivals  which  were  read  to  the  nearest  0.1  second,  which  showed  an  ISC  residual 


of  loss  than  S.O  seconds,  and  for  which  P-wave  station  corrections  were  available 
(Dziewonsld  and  Anderson,  1983).  Travel  time  residuals  were  computed  as  above 
relative  to  CNHC.  The  travel  time  data  from  faultless  show  fast  travel  times, 
with  the  fastest  travel  times  to  the  northeast.  These  travel  times  are  more 
negative  than  the  surface  receivers.  The  faultless  explosion  was  detonated  well 
below  the  valley  fill  in  saturated  tuff,  and  therefore  is  not  affected  by  the  sedi¬ 
ment  delays  as  are  the  surface  receivers.  The  relative  travel  time  advance  for 
faultless  suggest  the  possible  existence  of  high  velocity,  shallow  mantle 
material  beneath  the  Tertiary  volcanic  center.  Spence  (1974)  has  reported  a 
similar  observation  associated  with  the  Silent  Canyon  Caldera  at  the  Nevada  Test 
Site.  The  required  velocity  difference  beneath  the  cauldron  compared  to  more 
"typical**  Great  Basin  mantle  material,  is  within  the  range  predicted  by  Jordan 
(1979),  which  would  result  from  compositional  variation  in  the  upper  mantle 
associated  with  basaltic  differentiation. 


KKPKKErfCES 


Bache,  T.C.,  W.J.  Best,  R.R.  Blandford,  G.V.  Buiin,  D.G.  Harkrider,  E.J.Herrin,  A. 
Rymll,  and  M.J.  Shore  (1981).  A  technical  assessment  of  seismic  yield  esti¬ 
mation,  DARPA-  NMR-8 1-02,  Defense  Advance  Research  Projects  Agency, 
Arlington,  Virginia. 

Bache,  T.C.,  (1982).  Estimating  the  yield  of  underground  nuclear  explosions. 
Bull.  Sriama.  Soc.  Am,  72  Part  B,  p  13 1-188. 

Burdick,  L.J.,  and  D.V.  Helmberger,  (1978).  The  upper  mantle  P  velocity  struc¬ 
ture  of  the  western  United  States,  J.  Gaophys.  Ras.,  B3,  1899-1712. 

Chang,  A.C.,  and  D.H.  von  Seggern  (1980).  A  study  of  amplitude  anomaly  and  m* 
Mas  at  LASA  subarrays,  /.  Gaophys.  Bos.,  85,  4811-4828. 

Der,  Z.A.,  T.W.  McElfresh,  and  A.  O'Donnell  (1983).  An  investigation  of  the 
regional  variations  and  frequency  dependence  of  anelaatic  attenuation  in 
the  mantle  under  the  Uni  .ed  States  in  the  .5-4  Hs  band,  Submitted  to  the 
Gaophys.  J.R.  astr.  soc. 

Dstewonski,  A.  and  D.  Anderson,  (1983).  Travel  times  and  station  corrections  for 
p-waves  at  teleseismic  distance,  submitted  to  J.  Gaophys.  Bas. 


Eaton,  J.P.,  (1963).  Crustal  structure  from  San  Francisco,  California,  to  Eureka, 
Nevada,  from  seismic  refraction  measurements,  J.  Gaophys.  Ras.,  68, 
5789-5606. 

Everden,  J.  (1971).  Magnitude  versus  yield  of  explosions,  J.  Gaophys.  Res.,  75, 
1028-1032. 

Evemden,  J.F.  (1977).  Regional  bias  in  magnitude  versus  yield  measurements: 
Its  explanation  and  mode  of  evaluation,  written  communication. 

Green,  P.E.,  R.A  Frosch.  and  C.F.  Romany  (1965).  Principles  of  an  experimental 
large  aperture  seismic  array  (LASA),  Proc.  of  IEEE,  S3,  1821-1833. 

Herrin,  E.  (1968).  1968  seismological  tables  for  P  phases.  Bull  Saismo.  Soc. 
Am,  S3,  1193-1241. 

Hill,  D.t  and  L.  Paldser  (1986).  Crustal  structure  between  Nevada  Test  Site  and 
Boise,  Idaho,  in  77ie  Earth  Beneath  the  Continents:  Am.  Ceophys.  Union, 
Gaophys.  Mon.,  10.  391-419. 

Jordan.  T.  (1979).  Mineralogies,  densities  and  seismic  velocities  of  garnet  iherzo- 
lites  and  their  geophysical  implications,  in  The  Mantle  Sample:  Inclusions 
in  Kimberlites  and  Other  Volcanxcs,  edited  by  F.R.  Boyd  and  H.O.A.  Meyer, 
pp  1-14,  AGU,  Washington,  D.C. 

MacDonald,  GJL  (1969).  Composition  and  origin  of  Hawaiian  lavas,  in  studies  in 
vulcanology— a  memoir  in  honor  of  Howel  Williams:  GeoL  Soc.  America 
Mem,  116.  477-622. 

Pomeroy,  P.W.  (1980).  Seismic  discrimination  at  regional  distances,  Annual 
Technical  Report,  AFOSR  Contract  # F49620-80-C-0021 

Priestley,  K.F.,  and  J.N.  Brune  (1978).  Surface  waves  and  the  structure  of  the 
Great  Basin  of  Nevada  and  western  Utah,  J.  Gaophys.  Res.,  83,2285-2272. 

Priestley,  K.F.,  AS.  Ryall,  G.S.  Essie  (1982).  Crustal  and  upper  mantle  structure 
in  the  northwest  Basin  and  Range  Province,  Bull  Seismn.  Soc.  Am,  72, 
911-923. 

Scott,  D.H.  and  N.J.  Trash  (1971).  Geology  of  the  lunar  crater  volcanicfleld.  Nye 
County,  Nevada:  U.S.  Ceol.  Survey  Prof.  Paper  599-1,  22p. 

Spence,  W.  (1974).  P-wave  residual  differences  and  inferences  on  an  upper  man¬ 
tle  source  for  the  silent  oanyon  volcanic  centre,  southern  Great  Basin, 
Nevada,  Gaophys.  J.R.  a sir.  Soc.,  SB,  505-523. 

Springer,  D.L.  and  R.L.  Kinnaman  (1971).  Seismic  source  summary  for  U.S. 
under  ground  explosions,  1981-1970,  Bull  Saismo.  Soc.  Am,  81,  1073- 

1098. 


FIGURE  CAPTIONS 

Figure  1.  Location  map  of  Hot  Creek  Valley  in  the  vicinity  of  the  Faultless  shot- 
point  Stations  CNBC,  CNEC,  CNHM,  CNHR.  CNHC.  CNSB,  and  CNPS  denote  U.S. 
Geological  Survey  stations,  data  from  which  are  being  used  in  the  P -delay  study. 
The  average  relative  delay  is  given  following  the  station  designation.  Stations 
MOREY  MINE,  GROUND  ZERO,  and  NEEDLES  denote  University  of  Nevada  digital 
stations,  data  from  which  are  being  used  in  the  waveform  and  amplitude  studies. 
The  darkest  geologic  units  in  the  vicinity  of  CNEC  and  southward  denote  recent 
volcanic  deposits. 

figure  2.  Digital  seismograms  of  the  two  events  discussed  in  the  text. 

figure  3.  Comparison  of  the  power  spectra  for  vertical  component  seismograms 
showm  in  figure  2. 

Figure  4.  Left  -  Relative,  teleseismio  P-wave  residuals  as  a  function  of  azimuth; 
right  -  Relative  residuals  on  equal  area  lower  hemisphere  stereographic  projec¬ 
tions;  (  )+0.06-»+0.15  sec,  (  )+0.04-»— 0.04  sec,  (  *  )—  0.05-*— 0.15  sec. 
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ABSTRACT 

For  most  areas  in  the  western  Great  Basin  focal  mechanisms 
show  a  consistent  pattern  of  primarily  strike-slip  motion  for  shal¬ 
low  events  and  oblique  or  normal  slip  for  deeper  events.  However, 
orientation  of  the  axis  of  least  principal  stress  (T-axis)  is  different 
for  different  areas:  NW-SE  for  western  Nevada  and  the  Mono  Lake 
region,  and  NE-SW  for  the  Mammoth  Lakes  area.  Along  the 
remainder  of  the  Sierra  Nevada  frontal  fault  zone,  T-axes  show 
both  orientations. 

In  general  the  change  in  mechanism  with  depth  is  interpreted 
as  the  result  of  increasing  overburden  pressure,  resulting  in  rota¬ 
tion  of  the  maximum  compressive  stress  (P-axis)  from  horizontal 
at  depths  less  than  about  8  km  to  vertical  at  depths  greater  than 
about  9  km.  The  absence  of  normal-slip  events  at  depths  greater 
than  10  km  in  one  area  (Mono-Excelsior-Luning  zone)  may  be 
explained  by  a  larger  horizontal  compressive  stress  compared  with 
areas  that  do  have  normal  faulting  at  such  depths.  In  some  areas 
conjugate  right-  and  left-lateral  shear  on  nearly  vertical  fractures 
may  be  associated  with  the  formation  of  clusters  of  magma-filled 
dikes  at  shallow  depths. 

Assuming  strike-slip  faulting  to  be  characteristic  of  earth¬ 
quakes  with  depth  less  than  8  km  and  normal  faulting  for  events 
deeper  than  10  km,  extrapolation  of  measured  shear  stress  in  the 
upper  few  kilometers  of  the  crust  provides  a  rough  estimate  of  the 
maximum  and  minimum  principal  stress  as  a  function  of  depth.  At 
about  3  km  depth  we  find  that  St  and  S’ a  are  both  horizontal,  with 
values  of  about  104  and  55  MPa,  respectively.  At  depth  of  20  km  5  j 
is  vertical  and  equal  to  the  overburden  pressure,  about  530  MPa; 
S s  is  horizontal  and  about  280-300  MPa. 


INTRODUCTION 


The  Great  Basin  is  a  region  of  active  crustal  spreading  characterized  by  a 
system  of  horsts  and  grabens,  thin  crust,  high  heat  flow,  low  Pn  velocity  and 
relatively  low  gravity  (Stewart,  1971;  Thompson  and  Burke,  1974).  Structurally 
the  region  is  interpreted  in  terms  of  block-faulting,  with  alternating  grabens  and 
horsts  forming  a  series  of  parallel  valleys  and  mountain  ranges.  In  some  cases 
valleys  and  mountain  ranges  are  formed  respectively  by  the  downslope  and 
upslope  parts  of  tilted  blocks  (Stewart,  1978).  Tectonically  the  province  has 
been  formed  by  lithospheric  extension  in  an  almost  east-west  direction,  possibly 
involving  the  brittle  fragmentation  of  a  crustal  slab  above  a  plastically  extending 
lower  medium  (Stewart,  1971).  Different  tectonic  theories  to  explain  the 
regional  extension  are  given  by  Stewart  (1978). 

The  Great  Basin  is  characterized  by  a  complex  system  of  faults.  The  main 
fault  segments  have  north  to  northeast  trends;  northwest  trends  are  also 
observed,  primarily  in  the  southwest  part  of  the  region.  Both  normal  and 
strike-slip  motion  are  associated  with  Basin  and  Range  faulting.  The  Walker 
Lane,  a  northwest-trending  zone  in  western  Nevada,  is  characterized  by  right- 
lateral  strike  sUp  and  represents  a  boundary  between  NNW-trending  structures 
to  the  west  and  NNE  trends  to  the  east  (Slemmons  wt  aL,  1979;  Bell  and  Slem- 
mons,  1979).  On  the  other  hand  the  Olinghouse,  Midas,  Carson  and  Mono- 
Excelsior  lineaments  in  western  Nevada  have  northeast  trend  and  left-lateral 
strike  slip  (Sanders  and  Slemmons,  1979).  The  proportion  of  strike-slip  to  dip- 
slip  motion  appears  to  increase  from  north  to  south  in  the  region. 

Published  focal  mechanisms  for  the  western  Great  Basin  and  adjacent 
Sierra  Nevada  show  normal-,  oblique-  and  strike-slip  movement  (e.p.,Ryall  and 
Malone,  1971;  Koizumi  et  a L,  1973;  Smith  and  Lindh,  197B;  VanWormer  and  Ryall, 
1980).  Fault-plane  solutions  for  earthquakes  in  the  western  Great  Basin  indicate 
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a  generally  consistent  trend  (about  N60“W)  for  the  axis  at  minimum  compressive 
stress  ( T-asris ),  in  agreement  with  regional  lithospheric  extension  in  the  same 
direction  (Ryail  and  Malone,  1971).  However,  for  events  in  the  Sierra  Nevada 
frontal  fault  zone  (VanWormer  and  Ryall,  1980;  Cramer  and  Toppozada,  1980)  T- 
axis  trends  of  E-W  to  ENE-WSW  are  reported. 

Fault-plane  solutions  for  recent  earthquakes  in  northwestern  and  western 
Nevada  show  primarily  dip-  or  oblique-slip  movement,  and  surface  ruptures  show 
the  same  type  of  displacements.  The  Fairview  Peak  earthquake  in  1954,  for 
example,  had  a  focal  mechanism  with  a  strong  strike-slip  component,  agreeing 
very  well  with  geodetic  measurements  (Romney,  1957;  Slemmons,  1957).  The 
mechanism  for  the  Rainbow  Mountain-Stillwater  earthquakes  of  1954  also  had  a 
component  of  strike-slip  motion,  although  the  surface  ruptures  showed  pri¬ 
marily  vertical  slip  (Tocher.  1955;  1958).  Gumper  and  Scholz  (1971)  presented 
fault-plane  solutions  for  seven  large  earthquakes  plus  two  composite  solutions 
for  small  events,  all  located  in  the  western  Great  Basin  along  a  zone  from  the 
northern  Owens  Valley  to  Winnemucca.  They  found  predominantly  normal  fault¬ 
ing,  frequently  with  a  right-lateral  strike-slip  component.  They  interpreted  the 
latter  as  a  superposition  of  right-lateral  shear  stress  on  a  general  E-W  crustal 
extension.  As  an  exception,  the  authors  found  primarily  left-lateral  movement 
on  an  E-W  strike-slip  fault  for  a  composite  mechanism  in  the  Excelsior  Moun¬ 
tains,  and  interpreted  this  as  a  transform  fault  between  two  extensional  zones. 
However,  Ryall  and  Priestley  (1975)  reinterpreted  that  mechanism  as  oblique 
flip  on  a  fault  striking  NE. 

According  to  Ryall  and  Malone  (1971)  focal  mechanisms  in  the  Fairview 
Peak-Dixie  Valley  rupture  zone  are  consistent  with  an  interpretation  of  simple 
block  faulting,  with  faults  of  different  orientation  having  the  same  slip  direction 
and  a  consistent  orientation  of  about  N80°W  for  the  T-axis.  However,  Ryall  and 


-3- 


Vetter  (1982)  found  that  strike-slip  and  oblique-slip  mechanisms  tended  to 
occur  at  different  depths  in  the  Fairview  Peak-Dixie  Valley  zone,  and  attributed 
the  change  to  a  rotation  of  the  maximum  compressive  stress  ( P-axis )  related  to 
Increasing  overburden  pressure  with  depth.  Gumper  and  Scholz  (1971)  and 
VanWormer  and  Ryall  (1980)  found  that  mechanisms  south  of  the  Fairview  Peak 
zone,  in  Cedar  Valley,  Mina  and  Mono  Valley  were  consistent  with  right-oblique 
slip  on  faults  striking  NNW  and  dipping  east. 

In  studies  of  the  1980  Mammoth  Lakes,  California,  earthquake  sequence 
different  mechanisms  for  the  same  events  have  been  determined  from  local  and 
teleseismic  data.  For  the  larger  events  of  this  sequence,  fault-plane  solutions 
determined  from  recordings  in  the  California-Nevada  region  (Cramer  and  Toppo- 
zada.  1980;  Ryall  and  RyalL  1981a)  indicate  primarily  strike-slip  motion  (left- 
lateral  on  N-S  striking  planes  or  right-lateral  on  E-W  planes).  However,  long- 
period  data  collected  at  teleseismic  distance  for  the  same  events  require 
oblique  slip  on  moderately-dipping  faults  (Given  ef  at.,  1982)  or  a  non-double 
couple  component  in  the  moment  tensors  (Barker  and  Langston.  1982;  B.  Julian, 
written  communioation).  Given  et  at.  suggest  that  this  discrepancy  may  be  due 
either  to  distortion  of  the  radiation  pattern  by  structure  or  to  source  complex¬ 
ity.  In  the  latter  case,  they  hypothesize  that  the  short-period  arrivals  may 
reflect  failure  (strike-slip)  of  an  asperity  which,  when  broken,  allows  the  regional 
strain  to  be  relieved  with  the  mechanism  (oblique  slip)  determined  from  the 
long-period  data.  Observations  of  surface  ruptures  for  the  1980  Mammoth  Lakes 
earthquakes  indicate  primarily  dip  slip  or  oblique  slip  (Taylor  and  Bryant,  1980; 
Clark  et  a L,  1982).  Savage  and  Clark  (1982)  suggest  reinflation  of  a  magma 
chamber  at  mid-crustal  depth  in  Long  Valley  caldera  to  explain  both  observed 
surface  deformation  associated  with  the  1980  sequence,  as  well  as  the  earth¬ 
quakes  themselves. 
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A  detailed  discussion  of  the  oblique-slip  earthquakes  is  given  by  Zoback  and 
Zobaok  ( 1980a, b)  in  which  they  attribute  the  mechanism  to  slip  on  pre-existing 
faults  that  have  orientation  different  from  that  dictated  by  the  present  stress 
regime.  In  their  study  of  fault  patterns  in  north-central  Nevada  they 
hypothesize  a  stesss  field  in  the  form  Sy  =  Sna  >  Swtrw  (ie.  .Si  =  Sg  >  St),  so 
that  local  variations  in  the  magnitude  of  Sy  and  Sggg  determine  which  type  of 
faulting  occurs:  if  Sy  >  Sjour  the  slip  is  normal  or  oblique,  but  if  Sjukt  >  Sy  the 
movement  is  strike-slip.  Zoback  and  Zoback  also  suggest  that  the  Sierra  Nevada 
marks  a  transition  from  primarily  strike-slip  motion  along  the  San  Andreas  fault 
system  to  primarily  extensional  tectonics  of  the  Basin  and  Range  province.  In 
this  transition  zone  they  postulate  that  the  extensional  tectonics  and  high  heat 
flow  of  the  Great  Basin  extend  50-100  km  into  the  Sierra. 

A  number  of  authors  («.  g. ,  Wallace,  1979)  have  suggested  that  faulting  in 
the  Basin  and  Range  province  is  listric,  with  a  detachment  surface  at  15-20  km 
depth  (maximum  depth  of  earthquakes  in  this  region)  separating  the  brittle 
upper  crust  from  a  ductile  lower  crust.  However,  at  this  depth  movement  would 
take  place  through  very  weak  earthquakes  or  the  beginning  of  aseismic  creep, 
and  it  would  be  difficult  to  find  seismic  evdence  for  such  a  detachment  surface. 

In  a  previous  report  (Ryall  and  Vetter,  1982)  we  noted  that  focal  mechan¬ 
isms  for  earthquakes  in  the  Fairview  Peak-Dixie  Valley  rupture  zone  in  north- 
central  Nevada  were  consistent  with  NW-SE  extension,  as  reported  by  a  number 
of  previous  authors.  However,  there  appeared  to  be  a  correlation  between  the 
ratio  of  dip-slip  to  strike-slip  motion  and  Increasing  focal  depth,  such  that  shal¬ 
low  earthquakes  were  primarily  strike-slip  and  normal  faulting  predominated  at 
greater  depth.  This  is  similar  to  an  observation  by  Klein  si  aL  (1977)  that 
mechanisms  for  a  1972  Reykjanes  Peninsula  swarm  were  strike-slip  for  depth 
less  than  3  km  and  mostly  normal  for  greater  depth. 


In  the  present  investigation  we  have  determined  focal  mechanism  for  a 
large  number  of  earthquakes  in  the  western  Great  Basin  and  eastern  Sierra 
Nevada,  and  find  that  the  correlation  between  depth  and  mechanism  that  we 
observed  for  the  Fairview  Peak-Dixie  Valley  zone  generally  obtains  for  the  entire 
western  Great  Basin-eastern  Sierra  Nevada  region. 

ANALYSIS  AND  RESULTS 

The  data  base  for  this  study  consisted  mainly  of  recordings  from  the  Nevada 
seismic  network  of  30-50  stations  (depending  on  time  of  occurrence  of  a  particu¬ 
lar  event).  Supplementary  data  were  obtained  from  stations  operated  in  Califor¬ 
nia  by  the  University  of  California,  the  US  Geologioal  Survey  and  California  Insti¬ 
tute  of  Technology,  as  well  as  stations  in  Utah  operated  by  the  University  of 
Utah. 

Most  of  the  hypocenters  for  these  events  were  calculated  with  the  HYPO 71 
algorithm  of  Lee  and  Lahr  (1975);  a  few  were  determined  using  a  program  writ¬ 
ten  by  ¥.  A.  Peppin  for  Nevada  network  data.  In  both  cases  a  simple  crustal 
model  was  used  consisting  of  a  88  km- thick  layer  with  P-wave  velocity  8.0 
krn/sec  over  a  half-space  with  velocity  7.9  km/ sec.  A  north-trending  record  sec¬ 
tion  for  earthquakes  near  Mammoth  Lakes,  California,  indicated  that  a  better 
crustal  model  for  the  western  Great  Basin  would  include  an  intermediate  layer 
with  P-wave  velocity  8.7  km/sec  in  the  depth  range  22-28  km.  Work  by  other 
authors  (Thompson  and  Burke,  1974;  Prodehi,  1979;  Priestley  ef  a L.  1982)  indi¬ 
cates  that  the  crust  is  thicker  and  more  complex  along  the  Sierra  Nevada-Great 
Basin  boundary  zone.  However,  because  moat  of  our  observations  were  made 
either  in  the  Great  Basin,  where  intermediate  phases  do  not  appear  as  first 
arrivals  for  the  range  of  event  depths  in  this  study,  or  at  Pn  distances  in  Califor¬ 
nia,  the  simple  crustal  model  was  assumed  to  be  adequate  for  both  hypocenter 
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determination  and  analysis  of  first-motion  patterns. 

Most  of  the  events  used  in  this  investigation  had  location  quality  of  "A"  or 
**B",  usually  corresponding  to  standard  error  (ie..  precision  of  the  calculation) 
of  less  than  a  kilometer  in  depth  and  permitting  a  comparison  between  events  at 
different  depth  ranges  in  the  same  area.  However,  it  shoult  be  kept  in  mind  that 
absolute  aeouraoy  of  the  hypoeenter  calculation  is  a  function  of  the  model  used, 
and  the  depth  of  ail  events  in  a  given  area  could  vary  systematically  by  as  much 
as  a  few  kilometers. 


Depth  Distribution 

Figure  1  shows  frequency  of  earthquake  occurrence  as  a  function  of  depth 
for  all  A  and  B  quality  events  reoorded  and  located  since  1970.  The  plot  indi¬ 
cates  maximum  occurrence  in  the  7-10  km  range,  with  an  approximately  sym¬ 
metric  decrease  for  smaller  and  larger  depth.  Plotting  the  same  earthquakes  as 
a  function  of  both  depth  and  magnitude  (Figure  2)  shows  that  the  small  earth¬ 
quakes  scatter  over  a  greater  depth  range  than  the  stronger  events.  Thus, 
earthquakes  with  Jf&  2  2.0  are  strongly  concentrated  in  the  5-10  km  depth 
range,  while  moat  events  with  Mi  <  2.0  occur  at  depths  of  5-15  1cm  and  are  con¬ 
centrated  in  the  10-15  km  range. 

Focal  Mechanisms 

More  than  130  fault-plane  solutions  were  determined  for  this  study,  and  108 
were  sufficiently  well-determined  for  a  detailed  comparison  of  mechanism  as  a 
function  of  event  location  and  depth.  All  of  the  solutions  were  for  individual 
events  and  are  based  on  at  least  15  first-motion  observations;  in  many  cases  20- 
30  points  were  available.  Nearly  all  of  the  mechanisms  are  consistent  with 
extansional  (normal  or  oblique  slip)  or  strike-slip  faulting,  primarily  on  planes 
striking  in  a  northerly  direction  but  ranging  from  NW-SE  to  NE-STT.  In  agree- 
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msnt  with  our  previous  observations  for  the  Fairview  Peak-Dude  Valley  zone 
(Ryail  and  Vetter,  1982)  we  find  indications  for  a  systematic  change  in  focal 
mechanism  with  depth  for  nearly  all  parts  of  the  region  studied.  The  general 
picture  shows  strike-  and  oblique-slip  faulting  to  be  characteristic  in  the  upper¬ 
most  8-10  km  in  the  crust,  and  oblique  and  normal  faulting  to  predominate 
below  this  depth.  In  the  following  sections  mechanisms  will  be  discussed  for  four 
separate  areas  within  the  western  Great  Basin:  Mammoth  Lakes,  western 
Nevada,  the  Mono-Excelsior-Luning  zone,  and  the  Sierra  Nevada-Great  Basin 
boundary  zone  exclusive  of  the  Mammoth  Lakes  area. 

Mammoth  Lakma  area.  The  most  significant  activity  in  terms  of  number 
and  size  of  events  has  occurred  in  the  Mammoth  Lakes  area  since  1978;  in  gen¬ 
eral  the  Mammoth  Lakes  earthquakes  also  had  good  station  coverage  and  azimu¬ 
thal  distribution  of  observations.  We  studied  44  earthquakes  in  this  area,  36  of 
which  had  A  or  B  quality  hypocenter  determinations.  In  Figures  3  and  4  we  show 
fault-plane  solutions  for  34  of  these  events,  separated  respectively  into  groups 
with  focal  depth  less  than  or  greater  than  9  km.  All  the  plots  on  these  and  suc¬ 
cessive  figures  are  lower-hemisphere,  equal-angle  projections,  with  compression 
quadrants  shaded.  A  list  of  the  locations  and  mechanism  parameters  for  these 
events  is  given  in  Table  1.  The  earthquakes  in  the  figures  are  numbered  with 
increasing  depth,  in  accordance  with  the  listing  in  Table  1. 

Figure  3  indicates  that  for  Mammoth  Lakes  events  with  depth  less  than 
about  9  km  focal  mechanisms  are  rather  consistent  and  predominantly  strike- 
slip.  The  picture  is  less  clear  for  greater  depth  (Figure  4),  but  most  solutions 
show  a  relatively  larger  component  of  normal  slip  and  some  events  show  almost 
i  pure  normal  faulting.  There  are  some  strike-slip  events  at  depth  greater  than  9 

I 

km,  but  almost  no  shallow  shocks  were  found  with  predominant  dip-slip  motion. 
One  event  at  shallow  depth  even  had  a  reverse-faulting  component  (#26).  The 
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change  in  mechanism  with  depth  is  illustrated  by  Figure  5,  which  shows  two 
earthquakes  that  occurred  at  the  same  epicenter,  nine  minutes  apart  in  time 
and  with  a  difference  of  6  km  in  depth.  The  first,  with  a  strike-slip  mechanism, 
was  the  shallower  event  (depth  8.2  km);  the  second,  deeper  (14.2  km)  shock,  had 
a  strong  normal-slip  component. 

Western  Nmtada.  Figures  6  and  7  show  fault-plane  solutions  in  western 
Nevada,  respectively  for  depths  less  than  and  greater  than  9  km;  hypocenter 
coordinates  and  parameters  of  the  fault-plane  solutions  for  these  earthquakes 
are  listed  in  Table  2.  Events  in  the  Mono  basin.  Excelsior  Mountains  and  Luning 
area  will  be  shown  in  another  figure  and  discussed  separately.  Because  the 
western  Nevada  sample  contained  only  events  with  smaller  magnitude  than  the 
Mammoth  Lakes  sample,  station  coverage  for  the  former  was  not  as  complete  as 
for  the  latter,  and  focal  depths  were  not  as  well  constrained.  Only  the  best  solu¬ 
tions  are  shown  on  Figures  6  and  7,  but  many  of  them  had  C-quality  locations.  In 
general,  focal  depths  in  this  part  of  the  western  Great  Basin  are  somewhat 
larger  than  those  near  Mammoth  Lakes,  most  having  depth  greater  than  9  km. 
Similar  to  Mammoth  Lakes,  there  is  a  tendency  for  the  larger  shocks  to  cluster 
in  a  small  depth  range:  Moderate  earthquakes  In  the  Fairview  Peak  area  were 
concentrated  in  the  9-14  km  depth  range,  and  the  deepest  shocks  (depths 
around  20  1cm)  were  located  near  Schurz,  north  of  Walker  Lake.  Most  of  these 
deep  events  were  too  small  for  detailed  fault-plane  solutions;  of  the  few  events 
that  were  well-enough  constrained  (Figure  7)  all  had  oblique  or  normal  slip. 

In  spite  of  somewhat  less  constraint  on  focal  depth,  a  comparison  of  fault- 
plane  solutions  for  western  Nevada  gives  similar  results  as  for  the  Mammoth 
Lakes  events.  The  mean  depth  for  strike-slip  events  is  7  km,  while  for  oblique- 
and  normal-slip  events  it  is  12  and  13  km,  respectively. 

Mono-Exctlnor-Luning  zone.  The  most  active  part  of  the  Mono  basin- 
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Excelsior  Mountains-Luning  zone  La  just  east  of  Mono  Lake  near  the  California' 
Nevada  border.  Migration  of  earthquakes  along  this  zone  in  I960  was  inter¬ 
preted  by  C.  tide  (personal  communication)  as  a  propagating  strain  pulse,  sug¬ 
gesting  structural  continuity  at  depth.  As  a  result  events  in  this  zone  are  shown 
separately  on  Figure  8  for  the  Mono  basin  and  on  Figure  9  for  the  Excelsior 
Mountains  and  Lunin g  area.  Locations  and  parameters  of  the  fault-plane  solu¬ 
tions  are  listed  for  these  events  in  Table  3;  most  of  the  hypocenter  determina¬ 
tions  were  of  quality  A  or  B. 

Within  this  zone  the  tendency  for  normal  faulting  is  weak.  There  were  no 
normal-slip  events  at  all  in  the  Mono  basin  (Figure  8)  and  southern  Excelsior 
Mountains  (area  SW  of  Mina  on  Figure  9),  despite  the  fact  that  more  than  a  third 
of  the  events  had  depth  in  the  9-12  km  range.  There  is.  however,  a  trend  to  nor¬ 
mal  faulting  for  shoaks  in  the  northern  Excelsior  Mountains  and  Luning  area 
(areas  W  and  N  of  Mina  on  Figure  9).  In  agreement  with  other  areas,  the  mean 
depth  of  strike-slip  events  is  about  5  km  and  for  oblique-slip  events  it  is  9.5  km. 

Sierra  Nsvada-CrstU  Basin  boundary  sons.  Earthquakes  for  the  Sierra 
Nevada-Great  Basin  boundary  zone  (SNGBZ)  exclusive  of  the  Mammoth  Lakes 
area  are  shown  on  Figure  10.  Strictly  speaking,  only  those  earthquakes  on  the 
left  side  of  the  figure  are  in  the  Sierran  frontal  fauit  zone;  the  others  are  in  the 
westernmost  Great  Basin  but  have  characteristics  similar  to  events  along  the 
Sierra.  Two  of  the  mechanisms  (#4,  #8)  were  also  shown  on  Figures  6  and  7. 

Of  the  15  mechanisms  shown  on  Figure  10  and  listed  in  Table  4,  only  two  are 
tor  events  shallower  than  9  km.  In  general  a  tendency  toward  normal  faulting 
dominates  this  picture,  which  is  in  agreement  with  the  greater  depth  of  the 
earthquakes.  However,  focal  mechanisms  for  this  zone  are  not  as  consistent  as 
in  the  other  regions  shown  on  previous  figures:  T-axes  for  the  SNGBZ  show  both 
the  NW-SE  trend  of  the  western  Nevada  region  and  the  NE-SW  trend  of  the  Mam- 
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moth  sequence. 

DISCUSSION 

Two  main  observations  result  from  this  investigation.  First,  moderate-to- 
strong  earthquakes  in  the  western  Great  Basin  and  eastern  Sierra  Nevada  tend 
to  cluster  in  the  5-10  km  depth  range.  Second,  focal  mechanisms  for  events  in 
this  region  are  characterized  by  a  change  from  strike-  to  oblique-  to  normal-slip 
with  depth. 

The  tendence  for  larger  events  to  cluster  in  a  limited  depth  range  may  be 
explained  by  recent  papers  of  Sibson  (1982)  and  Meissner  and  Strehlau  (1982), 
dealing  with  the  maximum  depth  of  intraplate  earthquakes  as  a  function  of  the 
transition  from  brittle  to  ductile  behavior  of  rocks,  which  in  turn  is  a  function  of 
the  temperature  regime  of  the  region.  Based  on  results  of  laboratory  experi¬ 
ments  on  quartz-bearing  rocks  with  water  content  ranging  from  dry  to  saturated 
they  concluded  that  maximum  shear  stress  will  occur  at  a  characteristic  depth 
for  a  particular  region.  Comparison  of  frequency-depth  plots  for  well-recorded 
Intraplate  earthquakes  gave  good  agreement  with  the  calculations:  the  most  fre¬ 
quent  and  largest  earthquakes  tended  to  occur  at  a  depth  corresponding  to  the 
calculated  maximum  stress.  In  our  oase  the  concentration  of  earthquakes 
between  8  and  12  km  depth  is  in  accordance  with  predicted  maximum  shear 
stress  for  an  extensional  regime,  wet  quartz  rheology  (water-saturated),  and 
heat-flow  between  1.7  and  2.0  HFU  (Meissner  and  Strehlau,  1982).  The 
occurrence  of  only  small  <  2)  earthquakes  at  depth  greater  than  about  12 
km  suggests  that  shear  resistance  at  such  depths  is  weak  in  this  region. 

1b  the  previous  sections  we  have  shown  that  much  of  the  western  Great 
Basin  is  characterized  by  a  systematic  change  in  focal  mechanisms  —  from 
strike-slip  for  shallow  earthquakes  to  oblique-  or  normal-slip  at  depths  greater 
than  about  9-10  km.  As  illustrated  by  Figure  11,  this  appears  to  be  due  pri- 


manly  to  a  change  in  plunge  of  the  P-axis.  Thus,  for  the  Mammoth  Lakes 
sequence  Figure  11(a)  shows  that  T-axes  are  close  to  horizontal  and  have  rather 
consistent  NE-SW  trend.  On  the  other  hand,  while  the  trend  of  the  P-axes  is  also 
constant,  the  plunge  changes  —  from  almost  horizontal  for  shallow  earthquakes 
to  almost  vertical  at  larger  depth.  Similar  relationships  can  be  seen  on  Figures 
11(b)  and  11(c)  for  the  western  Great  Basin  and  the  SNGBZ,  respectively. 

This  observed  change  in  mechanism  with  depth  appears  to  occur  gradually. 
The  mean  depth  for  all  of  the  purely  strike-slip  earthquakes  is  about  5  km,  for 
oblique-  but  mostly  strike-slip  motion  it  is  6.2  km.  for  events  with  about  equal 
amounts  of  horizontal  and  vertical  movement  it  is  10.6  km,  and  for  earthquakes 
with  mostly  normal  slip  it  is  13.6  km.  While  there  are  a  few  mechanisms  on  Fig¬ 
ures  3-10  that  do  not  fit  this  pattern,  the  fact  that  average  behavior  with  depth 
ia  systematic  for  a  number  of  regions  strongly  suggests  that  the  type  of  faulting 
for  a  given  earthquake  is  determined  by  regional  extension  plus  a  depth- 
dependent  stress  component,  and  not  by  pre-existing  faults  (Zoback  and 
Zobaok,  1980a,  b)  or  simple  block  faulting  (Ryall  and  Malone.  1971). 

X*Hma£ion  of  Principal  Stresses.  These  observations  can  be  used  to  esti¬ 
mate  the  magnitude  of  the  principal  stresses  with  depth.  For  pure  strike-slip 
movement  the  maximum  compressive  stress  is  horizontal  and  the  intermediate 
stress  is  vertical,  while  for  normal  faulting  this  situation  is  reversed.  Our  results 
indicate  that  the  first  stress  pattern  exists  at  depths  less  than  5-8  km,  and  that 
the  second  pattern  dominates  at  depths  greater  than  9-10  km.  For  intermedi¬ 
ate  depths  oblique  mechanisms  indicate  a  gradual  rotation  of  the  P-axis  f  om 
horizontal  to  vertical,  an  exchange  of  the  maximum  and  intermediate  stress 
axes.  The  stress  field  suggested  by  Zoback  and  Zoback  (1980a,  b),  with 
5 1  w  St  »  St  and  St  vertical,  seems  to  be  characteristic  for  events  with  depth 
greater  than  9  km,  where  we  find  normal  or  normal-oblique  solutions.  At  depths 
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leas  than  about  6  1cm  the  predominance  of  strike-slip  solutions  indicates  that 
S i  >  Sz>  S a.  with  Sz  vertical.  For  most  of  the  solutions  S$  is  close  to  horizon¬ 
tal;  however,  the  average  trend  of  this  axis  is  NE-SW  for  events  in  the  Mammoth 
Lakes  area,  NW-SE  for  the  western  Great  Basin,  and  scattered  for  the  SNGBZ 
north  of  Mammoth  Lakes  (Figures  11a,  b  and  o,  respectively).  In  the  following 
paragraphs  we  consider  pore  pressure  P0,  overburden  pressure  Sr  and 
differential  stress  (Si  -  S3)  to  roughly  estimate  the  principal  stresses  Si  and 
S a.  Results  of  this  estimation  are  plotted  on  Figure  12. 

The  overburden  pressure  Sr  Is  determined  by  the  equation  Sr  =  pgz  where 
the  average  density  p  M  2.7pm /cm, "for  crustal  rocks.  Pore  pressure  Pq  is 
assumed  to  be  hydrostatic.  Numerous  authors  have  suggested  that  pore  pres¬ 
sure  is  hydrostatic  to  at  ieast  mid-crustal  depths  (s .g„  Brace  1971;  1972). 
Raleigh,  mt  aL  (1978)  found  pore  pressure  of  17  MPa  at  2  km  depth  in  Colorado,  a 
value  that  is  approximately  equal  to  hydrostatic.  As  shown  by  Zoback  and 
Zoback  (1980a,  b)  the  least  principle  stress  S  a  must  be  greater  than  the  pore 
pressure. 

Measurements  of  maximum  shear  stress  t —  *  (Si  —  Sa)/2  in  the  upper 
flew  kilometers  of  the  crust  have  been  interpreted  by  Haimson  (1977),  McGarr 
and  Gay  (1978),  Zoback  and  Roller  (1979)  and  Zoback  e<  at.  (1977).  These 
papers  indicated  that  shear  stress  increases  approximately  linearly  with  depth, 
and  Is  lower  in  soft  rock  than  in  hard  rock.  For  the  latter  McGarr  (1980)  gives 
tha  relationship 

Tan  a  5.B7  ♦  6.37  s ,  (l) 

where  rH  is  in  MPa  and  s  is  in  km.  An  estimate  of  the  shear  stress  at  10  km 

depth  in  north-central  Nevada  by  Zoback  and  Zoback  (1980a)  was  84-77  MPa, 

assuming  hydrostatic  pore  pressure.  This  is  in  good  agreement  with  equation 

(1),  which  gives  69.4  MPa  for  a  depth  of  10  km.  Figure  12  shows  the  differential 
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stress  (Si  -  Sa),  taken  as  twice  the  value  of  t™..  in  equation  (1). 

As  noted  above,  mechanisms  for  depths  larger  than  about  9  km  sire  normal 
or  oblique.  From  Jaeger  and  Cook  (1969)  frictional  sliding  on  preexisting  faults 
will  ocour  at  a  critical  ratio  of  the  effective  stresses, 

=  [(p2  +  l)1"*  +  m]8.  (2) 

where  fi  is  the  coefficient  of  friotion.  From  data  presented  by  Byerlee  (1978)  fj. 
is  about  0.66  at  a  depth  of  10  1cm  and  0.63  at  20  km.  Using  the  relationship 
a  as  S  —  Pq,  between  effective  stress  <r,  total  stress  5  and  pore  pressure  P0, 
equation  (2)  leads  to 


«?  -  (S  i  -  P  o)  .  , 

S*  ~  C0*a  +  i)‘*+M]*  °*  (3) 

where  in  the  case  of  normal  faulting  S  t  =  Sy,  the  overburden  pressure.  Equa¬ 
tion  (3)  gives  Sa  =  147  MPa  at  depth  10  km  and  298  MPa  at  20  km.  These  values 
for  the  least  principle  stress  are  shown  as  points  (l)  and  (2)  on  Figure  12. 

Somewhat  lower  estimates  of  Sa  can  be  obtained  as  follows;  Noting  that  the 
average  depth  for  obUque  mechanisms  is  10.8  km,  we  set  the  overburden  pres¬ 
sure  Sy  equal  to  the  maximum  horizontal  stress  Sg  at  that  depth  (point  3  on 
Figure  12.  281  MPa).  Then  from  equation  (l)  with  5 1  =  Sy  we  can  estimate  5  a  at 
10.6  and  20  km,  respectively,  as  134  and  283  MPa  (points  4  and  5  on  Figure  12). 

For  depths  less  than  about  5  km  fault-plane  solutions  for  the  western  Great 
Basin  are  almost  pure  strike-slip.  For  this  type  of  faulting,  according  to  Sibson 
(1974)  the  overburden  pressure  Sy  will  be  the  intermediate  stress  5g  and 

Sr  -  s,  =  (4) 


If  we  assume  that  pure  strike-slip  faulting  occurs  at  3  km  depth,  then  equations 
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(1)  and  (4)  give  values  for  Si  and  S3,  respectively,  of  104  and  55  MPa.  In  this 
case  Si  is  the  maximum  horizontal  stress  Sg  and  a  line  connecting  points  (3) 
and  (7)  on  Figure  12  gives  a  rough  estimate  of  the  increase  in  this  stress  with 
depth. 

Variation  in  T-axis  trends.  As  shown  by  Figure  11,  T-axes  (axes  of 
minimum  compressive  stress  S3)  for  this  study  have  maximum  plunge  of  40°  but 
in  most  cases  are  almost  horizontal.  The  mean  value  of  the  T-axis  plunge  for  all 
of  the  events  studied  is  about  8®,  with  no  systematic  change  with  depth.  Trend 
of  the  T-axes  is  NW-SE  for  events  in  the  western  Great  Basin  and  Mono  basin,  with 
a  peak  in  the  direction  N50-70®W  (Figure  3a). 

For  the  Mammoth  Lakes  area  T-axes  trend  NE-SW,  with  only  5  of  44  events 
failing  to  show  this  trend  (Figure  3B).  Earthquakes  along  the  rest  of  the  Sierra 
Nevada  firontal  fault  system  show  a  range  of  trends  from  NW-SE  to  NE-SW,  indi¬ 
cating  that  this  zone  is  a  region  of  transition  between  two  stress  systems.  This 
pattern  generally  agrees  with  the  finding  of  VanWormer  and  Ryall  (1980)  that 
meohanisms  along  the  Sierran  frontal  fault  zone  are  different  than  those  in  cen¬ 
tral  Nevada.  Our  findings  also  agree  with  those  of  Zoback  and  Zoback  (1980b) 
for  the  northern  Sierra  Nevada,  but  not  for  the  southern  Sierra  where  they 
found  NW-SE  T-axes.  Hill  (1962)  suggested  that  the  NE-SW  T-axis  for  the  Mam¬ 
moth  Lakes  sequence  represents  an  anomaly  in  a  region  of  otherwise  NW-SE 
extension,  but  our  results  suggest  that  the  NE-SW  orientation  is  not  unusual  for 
this  zone. 

Variations  in  fault  dip.  A  decrease  in  dip  of  the  fault  plane  with  increas¬ 
ing  focal  depth  is  indicated  by  the  parameters  listed  in  Tables  1-4.  Thus,  for  the 
shallow  event  shown  on  Figure  5  both  the  fault  and  auxiliary  planes  have  a  dip  of 
80®  or  more,  while  the  two  planes  for  the  deeper  event  have  dips  of  64*  and  38*. 
Regardless  of  which  plane  is  taken  to  be  the  fault  plane,  this  would  signify  a 
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change  from  steep  dip  near  the  surface  tc  smaller  dip  at  larger  depth.  Simi¬ 
larly,  the  mean  dip  of  the  fault  and  auxiliary  planes  for  all  the  events  we  studied 
is  about  80“-90#  for  earthquakes  shallower  than  6  km,  70"-75o  for  depths  of  6-9 
km,  and  55<’-60a  for  events  deeper  than  9  km.  This  shows  a  decrease  in  fault- 
plane  dip  by  about  25s  from  the  surface  to  about  15  km  depth.  While  this  obser¬ 
vation  might  by  itself  be  taken  as  support  for  listric  faulting,  other  evidence 
indicates  that  this  is  not  the  case.  Nowhere,  for  example,  have  we  seen  the  very 
shallow  dips  that  would  be  associated  with  a  detachment  surface  at  midcrustal 
depth.  Further,  the  sense  of  motion  —  strike-slip  at  shallow  depth,  oblique  or 
normal  at  larger  depth  —  is  not  appropriate  for  listric  faulting. 

Lithospheric  extension  by  faulting  and  magma  injection.  While  our 
observations  are  not  consistent  with  listric  faulting,  they  do  support  the  idea 
that  lithospheric  extension  in  the  western  Great  Basin  involves  both  tectonic  and 
magmatic  processes.  To  illustrate  this  we  recall  that  most  of  the  deeper  earth¬ 
quakes  in  this  area  are  associated  with  oblique-to-normai  slip  along  faults  paral¬ 
lel  to  the  Sierra  Nevada  frontal  fault  system,  striking  N-S  to  NW-SE  (figure  4). 
However,  the  shallower  events  have  strike-slip  mechanisms  —  right-lateral  slip 
on  planes  striking  about  N75°W  and/or  left-lateral  slip  on  planes  striking  about 
N10°E  (Figure  3).  In  an  earlier  paper  Hill  (1977)  pointed  out  that  shallow  strike- 
slip  events  are  commonly  observed  in  volcanic  regions,  and  proposed  a  model  to 
explain  this  pattern  in  which  conjugate  shear  failures  on  nearly  vertical  frac¬ 
tures  accompany  the  formation  of  magma-filled  dikes.  Comparison  of  his  model 
with  focal  mechanisms  for  the  Mammoth  Lakes  sequence  suggests  that  strike- 
slip  events  in  and  a round  Long  Valley  caldera  may  be  associated  with  the  for¬ 
mation  of  clusters  of  NW-trending  fissures  or  magma-filled  dikes  (Figure  13; 
Ryall  and  Ryall,  1983).  At  greater  depth,  extension  would  take  place  mainly 
through  faulting  instead  of  dike  formation.  The  existence  of  magma  at  shallow 
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depth,  at  least  within  the  caldera,  is  indicated  by  teleseismic  P-delays  (Steeples 
and  Iyer,  1976),  possible  reflections  on  an  explosion  profile  (Hill,  1976),  lack  of 
S-waves  for  paths  through  the  caldera  (Ryall  and  Ryall,  1981b),  uplift  of  the 
resurgent  dome  (Savage  and  Clark,  1982)  and  spasmodic  tremor  in  a  small  area 
near  Mammoth  Lakes  (Ryall  and  Ryall,  1983). 

Such  a  tectonic/volcanic  process  would  be  consistent  with  an  earlier 
suggestion  by  Lachenbruch  and  Sass  (1978)  that  lithospheric  extension  in  the 
Basin  and  Range  province  takes  place  through  a  combination  of  normal  faulting 
and  magmatic  intrusion  of  the  brittle  crust.  According  to  those  authors,  vol¬ 
canic  centers  like  Long  Valley  caldera  exist  "because  they  are  at  places  where 
the  lithosphere  is  pulling  apart  rapidly,  drawing  up  basalt  from  below  to  fill  the 
void."  Our  observation  that  a  systematic  change  in  focal  mechanism  with  depth 
is  also  characteristic  of  other  parts  of  the  western  Great  Basin  suggests  that  a 
combination  of  faulting  and  magmatic  activity  may  be  fairly  common  for  this 
entire  region. 

Finally,  it  should  be  noted  that  our  analysis  was  based  primarily  for  earth¬ 
quakes  in  the  range  Mi  =  3-5,  and  does  not  necessarily  apply  either  to  major  (M 
>  7)  earthquakes  in  the  region  or  to  microearthquakes  with  Mi  <  3.  For  the 
former,  field  evidence  indicates  that  major  oblique-  or  normal-slip  events  at 
mid-crustal  depths  rupture  to  the  surface  to  produce  range-front  faults  with  the 
same  type  of  slip.  In  the  case  of  the  strong  earthquakes  at  Mammoth  Lakes  in 
1980,  it  is  possible  that  strike-slip  movement  associated  with  the  breaking  of  a 
shallow  asperity  may  have  triggered  a  larger  oblique-slip  event  at  depth,  as  sug¬ 
gested  by  Given  et  a l  (1982),  leading  to  oblique-slip  motion  observed  at  the  sur¬ 
face  by  Taylor  and  Bryant  (1980)  and  Clark  et  at.  (1982).  To  extend  this  study 
we  are  currently  studying  smaller  (Mi  <  3)  earthquakes  to  try  and  determine 
whether  they  behave  in  the  same  systematic  way  as  the  moderate-sized  events. 
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Table  2.  Western  Nevada  (Figures  6,  7).  Column  headings  same  as  Table 
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Table  3.  Mono  -  Excelsior  -  Luning  Zone  (Figures  8,  9).  Column  headings  same  as  Table 
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Table  4.  Sierra  Nevada  -  Great  Basin  Boundary  Zone  (Figure  10).  Column  headings  same  as 
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FIGURE  CAPTIONS 

Figure  1.  Earthquake  frequency  as  a  function  of  focal  depth,  for  events  with  A 
and  B  quality  hypocentral  determination  in  the  western  Great  Basin  and  eastern 
Sierra  Nevada,  1970-1981.  N  —  number  of  earthquakes. 

Figure  2.  Earthquake  frequency  (percent)  as  a  function  of  depth,  (a)  —  Magni¬ 
tude  ML  greater  than  2;  (b)  —  ML  less  than  2.  Only  A  and  B  quality  events  are 
shown. 

Figure  3.  Fault-plane  solutions  for  earthquakes  with  depth  less  than  9  km  in  the 
Mammoth  Lakes  area,  eastern  California.  (Lower-hemisphere,  equal-angle  pro¬ 
jection,  shaded  areas  are  compressions,  numbers  correspond  to  Table  l). 

Figure  4.  Fault-plane  solutions  for  earthquakes  with  depth  greater  than  9  km  in 
the  Mammoth  Lakes  area,  eastern  California.  Numbers  correspond  to  Table  1. 

Figure  5.  Fault-plane  solutions  for  two  Mammoth  Lakes  events  that  occured  on 
Aug.  9,  1981  at  the  same  epicentral  location.  Time  (GMT),  coordinates,  depths 
(D),  magnitudes  (M),  and  quality  (B)  are  given.  Closed  circles  —  compression, 
open  circles  —  dilatation,  T  —  T  axis,  P  —  P  axis. 

Figure  8.  Fault-plane  solutions  for  earthquakes  with  depth  less  than  9  km  in  the 
western  Great  Basin.  Numbers  correspond  to  Table  2. 

Figure  7.  Fault-plane  solutions  for  earthquakes  with  depth  greater  than  9  km  in 
the  western  Great  Basin.  Numbers  correspond  to  Table  2. 

Figure  8.  Fault-plane  solutions  for  earthquakes  east  of  Mono  Lake  at  the 
California-Nevada  border.  Numbers  correspond  to  Table  3. 
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Figure  9.  Fault-plane  solutions  for  earthquakes  in  the  Excelsior  Mountains  and 
Luning  area  of  western  Nevada.  Numbers  correspond  to  Table  3. 

Figure  10.  Fault-plane  solutions  for  earthquakes  in  the  Sierra  Nevada-Great 
Basin  boundary  zone.  Numbers  correspond  to  Table  4. 

Figure  11.  Lower-hemisphere  equal-angle  projection  showing  P-axes  (solid  cir¬ 
cles)  and  T-axes  (open  circles),  (a)  --  Mammoth  Lakes  area;  (b)  —  western 
Nevada  and  Mono-  Excelsior-Luning  zone;  (c)  —  Sierra  Nevada-Great  Basin  boun¬ 
dary  zone. 

Figure  12.  Estimated  principal  stresses  as  a  function  of  depth.  P 0  —  hydrostatic 
pore  pressure,  Sj  —  Sa  --  twice  the  maximum  shear  stress  measured  for  hard 
rocks  and  extrapolated  downwards,  S*=S 3  —  smallest  horizontal  principal 
stress,  Sy  —  overburden  pressure,  Sg  —  greatest  horizontal  principal  stress. 
Shaded  area  shows  range  of  S3.  Points  with  numbers  are  described  in  the  tex- . 

Figure  13.  Model  to  explain  strike-slip  faulting  for  shallow  Mammoth  lakes  earth¬ 
quakes.  Top  left  —  rose  diagram  showing  orientation  of  T  axes,  top  right  --  rose 
diagram  showing  strike  of  the  two  fault  planes,  bottom  —  model  with  vertical 
fissures  (dikes)  forming  together  with  conjugate  right-  and  left-lateral  shears  on 
planes  shown. 
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